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ABSTRACT 


The  determination  of  the  LaPlacian  of  a  mapped  quantity,  such  as 
the  geostrophic  vorticity,  is  an  important  but  a  time  consuming  process. 
The  common  procedure  involves  a  number  of  interpolations  and  an  algebraic 
computation  for  each  point  at  which  the  evaluation  is  made.  A  device 
is  presented  in  which  the  slow  process  of  interpolation  is  replaced  by 
the  faster  and  more  reliable  process  of  measurement  of  distance  between 
isopleths  of  the  mapped  quantity,  and  the  algebraic  calculation  is  re¬ 
placed  by  the  rapid,  automatic  operation  of  a  DC  analog  computer.  The 
numerical  result  is  presented  as  the  deflection  of  a  millivoltmeter. 

In  practice  the  use  of  this  computer  proves  to  be  a  rapid,  repro¬ 
ducible,  sufficiently  accurate  means  of  determining  the  geostrophic  vor¬ 
ticity  from  an  analysis  of  the  contours  of  an  isobar ic  surface. 


1.  Introduction 


The  operation  of  determining  the  two-dimensional  LaPlacian  of  a 
mapped  variable,  ^72h  *  9  2h/^  +  ^2h/dy^,  where  x  and  y  are  hori¬ 

zontal  Cartesian  coordinates,  is  a  recurring  problem  in  meteorology 
today.  As  the  application  of  theoretical  results  to  forecasting  grows, 
more  and  more  operational  weather  groups,  as  well  as  research  groups, 
find  that  they  must  devote  appreciable  effort  to  the  evaluation  of  this 
and  other  differential  expressions.  The  outstanding  example  is  the 
evaluation  of  the  geostrophic  vorticity  from  the  contours  of  an  analyzed 
isobaric  chart. 

The  usual  approach  to  the  numerical  evaluation  of  a  differential 
expression  is  to  approximate  derivatives,  which  are  ratios  of  infini¬ 
tesimal  increments,  with  ratios  of  finite  increments.  Thus  a  differ¬ 
ential  expression  is  reduced  to  an  algebraic  one.  For  example,  if  in 
fig.  1  the  values  of  the  variable  h  are  interpolated  at  the  five  points 
indicated, 

V2h  a  *  h2  *  ^  ■!) 

L2 

where  the  subscript  identifies  the  point  at  which  a  value  of  h  is  inter¬ 
polated,  and  L  is  the  constant  finite  distance  separating  points  1  to  ii 
from  point  5.  This  technique  requires  several  interpolations  and  an 
algebraic  computation  for  each  value  of  the  LaPlacian  which  is  calculated. 
Both  parts  of  the  procedure  are  time  consuming  and  subject  to  error; 
thus  the  process  is  not  completely  satisfactory. 

One  way  by  which  an  improvement  might  be  made  is  to  utilize  a  ratio 
of  finite  differences  in  which  the  increment  of  the  variable  h  is  kept 
fixed  and  equal  to  the  increment  between  successive  isopleths  in  a  map 
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Fiy.  1.  A  grid  of  lour  joints  equally  spacer  around  a  central  point. 
Interpolated  values  at  these  five  points  are  used  to  evaluate  the 
LaPlacian  of  the  field  of  h  by  means  of  equation  (1). 


Ah 


Fig.  2.  Distances  along  orthogonal  directions  separating  adjacent 
lines  of  the  variable  h.  These  distances  can  be  used  to  deter¬ 
mine  the  LaPlacian  of  the  field  of  h  by  means  of  equation  (3). 


3 


of  h.  The  distance  between  isopleths  would  then  vary,  and  distance 
measurement  would  replace  interpolation.  In  principle,  one  should  be 
able  to  measure  distances  more  rapidly  and  with  fewer  errors  than  one 
could  make  interpolations.  In  addition,  it  might  be  possible  to  gear 
the  mechanical  device  which  measures  the  distances  to  a  simple  computer 
which  would  do  the  algebra  necessary  to  evaluate  the  LaPlacian.  Thus 
both  of  the  cumbersome  steps  in  the  usual  process  could  be  replaced. 

In  subsequent  sections,  we  shall  show  one  way  in  which  these  suggestions 
can  be  carried  out. 

2.  A  finite  difference  approximation  to  the  LaPlacian 

In  fig.  2,  the  orthogonal  distances  separating  an  isopleth  of  h 
from  adjacent  ones  are  labelled  A,  B,  C,  and  D.  If  the  increment  in 
h  from  one  isopleth  to  the  next  is  the  constant  value  Ah,  approxi¬ 
mations  to  the  first  derivatives  in  the  x-direction  are  Ah/A  and 
Ah/C.  Correspondingly  the  approximations  to  the  first  derivatives 
in  the  y-direction  are  -  Ah/B  and  -  Ah/D.  The  second  derivatives 
are  approximately 

Ah  Ah  -  Ah  Ah 

A  “  C  B  +  D 

7" - -  and  — r“” - — 

i  (A  +  C)  i  (B  +  D; 

Thus, 

2*h  [(M)  rh*  (b-D  s-b] 

or 

V2h  «  2  Ah[I^?%j  +  Sbfa“;DD)  ] 


To  sirplify  this  expression,  we  multiply  and  divide  the  first  term  in 
the  brackets  by  A  +  C,  and  perform  the  same  operations  on  the  second  term 


h 


with  B  +  D.  The  result  is 


V^hiS  2Ah 


C2  -  A2  B2  -  D2 
AC  (A  ;►  C)2  +  BD(B  +  0)2 


(2) 


Now  a  quantity  like  (A  +  C)2  *  UAC  +  (A  -  C)2  may  be  approximated  quite 
well  if  A  and  C  are  not  greatly  different.  The  fractional  error  in 
assuming  that  (A  +  C)2  -  IiAC  is  (A  -  C)2/(A  +  C)2.  This  error  is  zero 
when  A  *  C  and  is  about  10  per  cent  when  A  ■  2C.  When  we  introduce 
this  kind  of  approximation  into  both  terms  on  the  right  of  equation  (2) 
and  simplify,  we  get 


V2h= 


Ah 

“T 


1  l  11 


(3) 


Equation  (3)  involves  the  sum  of  four  inverse  squares,  and  is  in  a 
convenient  form  for  the  design  of  an  electrical  analog  which  can  auto¬ 
matically  perform  the  algebra.  However,  it  is  not  the  only  approxima¬ 
tion  to  the  LaPlacian  which  might  prove  useful. 

When  equation  (3)  is  used  to  estimate  the  geostrophic  vorticity, 
the  approximations  mean  that  in  cases  of  marked  shear  the  magnitude  of 
the  vorticity  will  be  overestimated  somewhat. 

3.  Theory  of  an  electrical  analog 

If  equation  (3)  is  used  to  approximate  the  LaPlacian,  the  problem 
can  be  reduced  to  mechanical  measurement  of  the  four  distances  and 
electrical  calculation  of  a  sum  of  inverse  squares.  The  distances 
involved  may  be  measured  by  means  of  a  rack  controlled  by  a  gear.  Thus 
the  distances  are  proportional  to  rotation  of  the  gear.  The  gear  can 
be  attached  easily  to  one  or  more  variable  electrical  resistances.  Thus 
the  problem  is  to  find  a  suitable  measuring  circuit.  The  principle  of 
the  method  which  was  adopted  can  be  seen  by  considering  fig.  3.  The 


upper  portion  of  this  figure  shows  a  simple  voltage  divider,  in  which 
the  potential  across  the  fixed  resistance,  J^,  is  inversely  proportional 
to  the  first  power  of  the  variable  resistance  R.  If  one  thinks  of  ap¬ 
plying  another  similar  voltage  divider  as  in  the  bottom  part  of  fig.  3, 
the  output  will  be  approximately  proportional  to  an  inverse  square  of 
R.  Since  the  two  resistances,  R,  can  be  coupled  jointly  to  the  gear, 
it  is  possible  to  get  an  output  voltage  from  such  a  circuit  which  is 
inversely  proportional  to  the  square  of  the  distance  traversed  by  the 
rack.  This,  of  course,  is  the  desired  situation. 

The  double  voltage  divider  of  fig.  3  has  an  actual  output  voltage 
which  is  more  complex  than  the  approximation  given  above.  We  shall 
now  calculate  the  true  output,  £.  The  total  resistance,  R^,  of  the 
entire  circuit  is 


R  *  ♦  VCR  ♦  Rc) 


R2  +  3RRp  ♦  Rp2 
R  +  2Ro 


If  the  input  voltage  is  V,  the  total  current,  1^,  flowing  through  the 
circuit  is 


U  V  -  V&+2gttl.. 
1  Rt  r2  +  3RR0  +  V 


00 


Now  the  voltage  drop  from  points  A  to  B  in  fig.  3  is  IjRq  and  also 
I2(R  +  Rq).  Since  It  ■  I^  +  I2,  we  may  write  for  I2 

H  irrs^ 

Substitution  for  1^  from  (U)  yields 

I9  m  v  — * - ^0- - 

2  R^  +  3RRq  +  Rq2 

Thus  E  *  I2R0  way  be  written  as 

d  2 


.3.  At  the  topi  a  simple  voltage  divider  ill  which  the  output 
voltage  is  inversely  proportional  to  the  first  power  of  the 
variable  resistance,  R.  At  the  bottom  a  double  voltage  divider 
in  which  the  output  voltage  is  approximately  inversely  proportional 
to  the  square  of  the  variable  resistance,  R. 
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This  may  also  be  written  as 

E  -  V 


«o2 

(R  +  |  Ro)2  Ro2 


(5) 


If  the  last  term  in  the  denominator  of  (5)  can  be  kept  sufficiently  small 
compared  to  the  first  term,  we  may  write 

Ro2 

E«V - - - -  (6) 


in  which  the  output  voltage  is  approximately  proportional  to  the  inverse 
square  of  a  linear  function  of  R,  as  desired.  The  fractional  error  in 
this  approximation  is  5Ro^A(R  +  3Rq/2)2.  Expression  (6)  is  an  under¬ 
estimate  of  the  true  value  given  by  (5).  However,  there  is  little 
difficulty  in  keeping  this  error  at  an  acceptably  low  value.  Thus,  for 
the  adopted  value  R0  ■  1$  ohms,  if  one  does  not  allow  R  to  go  below 
1*2.5  ohms  the  error  does  not  exceed  5  per  cent.  The  restriction  to  R 
greater  than  about  1*0  ohms  is  readily  accomplished}  thus  the  maximum 
error  from  this  source  is  tolerable. 

The  simple  theory  given  above  is  for  one  measuring  circuit.  Four  of 
these  are  necessary,  connected  in  series  and  with  appropriate  polarities. 
The  polarities  are  such  that  the  errors  discussed  above  tend  to  cancel. 

1*.  Design  details 

The  completed  instrument  is  shown  in  use  in  fig.  1*.  The  four  racks 
are  attached  to  shafts  which  slide  in  grooves  machined  in  plexiglas.  The 
steel  racks  are  5  inches  long  and  are  operated  by  brass  gears  of  SO  pitch 
and  pitch  diameter  1.360  inches.  This  permits  a  range  in  the  displacement 
of  the  end  of  a  shaft  from  the  center  point  by  a  factor  of  10. 

Each  gear  is  coupled  to  a  wheel  for  manual  control  and  to  a  pair  of 
matched,  coupled,  variable  resistances  (R  in  fig.  3).  These  potentio¬ 
meters  have  a  maximum  resistance  of  500  ohms,  are  accurate  to  within  2 
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Fig*  U.  The  computer  in  operation.  At  the  lover  right  is  the  measuring  frame.  The 
central  hole  has  been  placed  on  a  contour  and  the  ends  of  the  measuring  shafts 
have  been  set  on  adjacent  contours.  The  variable  resistances,  gears,  and  con¬ 
trol  knobs  are  at  the  four  comers  of  the  frame*  At  the  upper  left  is  a  case 
containing  the  power  supply,  resistors,  etc*  On  the  face  of  the  case,  from 
left  to  right,  are  the  on-off  switch,  millivoltmeter,  and  latitude  selector 
switch. 
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per  cent  of  the  required  value  at  any  setting  between  hS  and  500  ohms, 
and  have  a  rotation  of  300*.  They  were  supplied  by  the  Clarostat  Co. 
Two-watt  potentiometers  were  ordered;  however,  it  would  have  been  con¬ 
venient  to  have  had  a  higher  wattage. 

The  plexiglas,  through  which  the  underlying  map  may  be  viewed, 
had  a  hole  drilled  to  mark  the  center  point.  Beside  this,  two  large 
openings  were  cut  to  permit  written  entries  to  be  made  on  the  map. 

The  remainder  of  the  circuitry,  the  power  supply  and  meter  are 
mounted  in  a  separate  case  connected  to  the  meter  by  a  twelve -conductor 
cable.  A  schematic  wiring  diagram  is  given  in  fig.  5.  The  15  ohm 
resistors  (Rg  of  fig.  3)  are  accurate  to  within  1  per  cent.  The  power 
supply  for  each  measuring  circuit  consists  of  three  mercury  batteries 
in  series,  giving  a  no-load  potential  of  ii.035  volts.  These  are  Mallory 
KM-1;2R.  cells  with  a  rated  capacity  of  lii,000  mi  Hi  ampere-hours.  Mercury 
cells  were  chosen  because  of  their  long  shelf  life  and  constancy  of  out¬ 
put.  In  practice,  under  load,  the  output  is  closer  to  3.75  volts.  In 
order  to  prevent  excessive  current  drain,  each  circuit  is  normal*,  open. 

To  make  a  reading,  a  switch  must  be  depressed  which  closes  all  four 
circuits  simultaneously. 

The  meter  used  to  measure  output  is  a  200-0-200  millivoltmeter  of 
5U7  ohms  resistance.  In  order  to  prevent  overshooting  of  the  pointer 
at  nearly  full  scale  deflections,  a  500  microfarad  capacitor  was  con¬ 
nected  in  parallel  across  the  meter.  Since  geostrophic  vorticity  for 
an  isobar ic  surface  is  %  *  ^  where  g  is  the  acceleration  of  gravity, 

f  is  the  Coriolis  parameter,  and  h  is  geopotential  height,  it  is  neces¬ 
sary  to  allow  for  the  variation  with  latitude  of  the  factor  l/f.  This 
is  done  by  providing  a  number  of  resistors  which  may  be  put  in  series 
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with  the  meter  to  cut  down  its  sensitivity  by  appropriate  amounts.  A 
selector  switch  inserts  the  resistor  appropriate  to  each  latitude  in  5* 
steps,  from  25*  to  65*. 

The  numerical  values  which  have  been  given  for  the  various  components 
represent  only  one  suitable  combination.  There  are,  of  course,  many 
more.  In  any  given  subsequent  case  these  values  would  have  to  be  re¬ 
calculated  to  comprise  a  consistent  set.  The  chief  requirements  are 
that  R0  be  kept  as  low  as  the  wattage  rating  of  the  potentiometers  R 
will  allow,  and  that  the  resistance  of  the  meter  be  as  high  as  is  con¬ 
sistent  with  convenience  and  sensitivity  needs. 

5.  Operational  procedure 

The  following  is  a  step-by-step  procedure  for  utilizing  this  com¬ 
puter  to  determine  the  geostrophic  vorticity. 

A.  The  measuring  frame  is  placed  on  the  map,  with  the  central  hole 

resting  on  a  contour  in  the  region  where  a  measurement  is  desired,  as  in 
fig.  1*.  The  frame  should  be  oriented  so  the  measuring  arras  make  an  angle 

of  about  1*5*  with  the  contours.  This  angle  is  not  critical.  In  tne 

model  which  has  been  built,  the  plug  on  the  measuring  frame  is  always 
kept  towards  low  heights. 

B.  The  latitude  control  switch  is  turned  to  the  latitude  position 
nearest  to  the  actual  value. 

C.  The  four  measuring  arms  are  adjusted  until  their  ends  rest  on 
the  contours  adjacent  to  the  central  one. 

D.  The  on-off  switch  is  depressed,  causing  the  meter  to  indicate. 

E.  The  meter  reading  is  recorded  on  the  map  through  one  of  the 

slots  in  the  plexiglas  provided  for  this  purpose. 
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In  practice,  two  situations  have  been  encountered  in  which  diffi¬ 
culty  is  experienced.  First,  if  the  contour  curvature  along  a  trough 
or  ridge  line  is  very  great,  the  measuring  arms  may  not  intersect  an 
adjacent  contour.  To  keep  this  trouble  to  a  minimum,  it  is  advisable 
to  use  a  base  map  of  sufficiently  small  scale.  When  the  difficulty 
arises,  it  may  be  reduced  further  by  changing  the  U5*  orientation  of 
the  measuring  arms  to  the  contours,  or  by  using  intermediate  contours. 
Whenever  intermediate  contours  (half  the  usual  height  increment)  are 
used,  the  meter  results  must  be  divided  by  two. 

The  second  difficulty  arises  in  closed  centers.  The  following  ap¬ 
proximate  procedure  has  been  used  in  such  situations  with  success. 

Place  the  central  point  of  the  measuring  frame  in  the  center  of  the 
closed  contour.  Put  two  of  the  measuring  arras  (having  the  same  polarity) 
on  the  innermost  closed  contour  with  the  other  two  arms  as  far  out  as 
possible.  Note  the  meter  reading,  then  run  out  the  two  arms  which  were 
in  and  place  the  two  arms  which  were  out  on  the  closed  contour.  Note 
the  meter  reading.  Next,  multiply  the  absolute  sum  of  the  two  readings 
by  that  fraction  of  the  standard  height  interval  between  contours  which 
it  is  estimated  lies  between  the  closed  contour  and  the  center  of  the 
system.  This  procedure  is  based  upon  the  fact  that  a  measuring  circuit 
gives  a  negligible  output  when  its  measuring  arm  is  as  far  out  as  possi¬ 
ble. 

In  certain  circumstances,  such  as  with  small  height  gradients,  it 
may  be  desirable  to  utilize  intermediate  contours,  in  which  case  the 
meter  reading  should  be  halved,  as  indicated  above.  When  the  gradients 
are  very  large,  it  is  advisable  to  utilize  every  other  contour  (doubled 
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interval).  Then  the  meter  readings  must  be  doubled  before  being  plotted. 
6,  A  comparative  test 

In  order  to  test  the  accuracy  and  utility  of  the  computer,  a  syn¬ 
thetic  height  field  was  constructed  as  shown  in  fig.  6.  This  is  the  stun 
of  two  fields,  one  of  circular,  equally  spaced  contours,  and  the  other 
of  straight  parallel  contours  extending  from  left  to  right  with  anti- 
cyclonic  shear.  The  true  LaPlacian  of  this  field  is  known  and  is  shown 
in  fig.  7a.  The  value  of  vorticity  at  the  center  of  the  pattern  is  in- 
finitej  thus,  it  is  inevitable  that  any  finite  difference  method  for 
estimating  the  LaPlacian  will  be  in  error  close  to  the  center.  However, 
it  will  be  of  interest  to  note  how  close  various  methods  can  come  to 
this  unattainable  pattern. 

The  LaPlacian  evaluated  by  means  of  the  computer  is  shown  in  fig. 

7b.  It  can  be  seen  that  the  general  shape  of  the  true  pattern  has  been 
portrayed  successfully.  The  smoothing  involved  in  the  use  of  finite 
differences  causes  the  absolute  values  to  be  smaller  than  they  should  be 
and  weakens  the  gradients,  especially  near  the  center.  Fig.  7c  gives 
the  LaPlacian  obtained  by  the  graphical  technique  of  evaluating  equa¬ 
tion  (1),  given  by  Fj^rtoft  (1952).  The  differentiation  distance,  L, 
was  taken  to  be  U*  of  latitude.  It  is  obvicus  that  this  method  smooths 
the  field  considerably  and  so  markedly  lowers  the  magnitudes  and  gra¬ 
dients.  Fig.  7c  could  hardly  be  accepted  as  an  adequate  representation 
of  fig.  7a,  except  for  the  most  gross  applications.  The  difference 
between  7b  and  7c  is  largely  due  to  the  change  in  the  distance  over 
which  the  finite  differences  are  taken  in  the  two  cases.  In  7b  the 
differentiation  distance  is  roughly  the  distance  between  the  contours 
of  fig.  6.  This  is  variable,  but  is  approximately  1*  of  latitude. 

It  is  desirable  to  compare  the  results  obtained  with  the  computer  to 


.  6.  Synthetic  contour  pattern  for  an  isobaric  surface.  The  lines 
are  labelled  with  values  in  tens  of  feet.  The  distance  oetween 
adjacent  tick  marks  represents  ii*  of  latitude. 


Fig.  7.  Four  different  representations  of  the  LaPlacian  of  the  height 
field  of  fig.  6.  a  is  the  true  LaPlacian  in  which  isopleths  have 
not  been  drawn  in  the  hatched  region)  3  is  the  field  determined  by 
the  graphical  method  of  Fj^rtoft  using  a  differentiation  distance 
of  U*  of  latitude)  D  is  the  field  determined  by  interpolation  using 
a  differentiation  distance  of  1*  of  latitude.  In  all  four  cases  the 
numerical  values  are  expressed  in  feet  and  are  directly  comparable. 

•\  The  values  are  those  of  ^h  where  L  is  of  latitude  in  each 

case.  In  B,  C>  and  D  the  lines  have  been  smoothed  to  a  roughly  equal, 
small  degree. 
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those  obtained  with  a  standard  technique  using  a  comparable  differenti¬ 
ation  distance.  Fig.  7d  gives  the  results  of  applying  equation  (1), 
where  L  *  1*  of  latitude.  Although  the  central  values  and  gradients 
are  satisfactory,  it  is  clear  that  the  field  of  the  LaPlacian  is 
highly  irregular  and  erratic.  This  is  because  equation  (1)  requires 
that  two  successive  differences  of  interpolated  values  be  taken.  As 
the  differentiation  distance,  L,  becomes  smaller  it  is  necessary  that 
the  accuracy  of  interpolation  increase,  in  order  that  the  accuracy  of 
the  differences  remain  constant.  When  the  grid  distance  is  as  small 
as  1*  of  latitude,  the  results  are  affected  significantly  by  random 
errors.  This  demonstrates  that  information  concerning  the  derivatives 
of  a  field  may  be  extracted  more  accurately  by  measuring  distances 
between  lines  in  the  field  than  by  interpolation. 

It  is  perhaps  worth  pointing  out  that  in  obtaining  fig.  7d  the 
graphical  process  could  not  be  used.  This  is  because  the  graphical 
technique  requires  that  the  difference  be  found  between  a  smoothed  field, 
(h^  +  h2  +  h^  +  h^)/U,  and  the  unsmoothed  field.  If  the  distance,  L, 
is  too  small,  there  is  so  little  difference  between  the  smoothed  and 
unsmoothed  fields  that  when  superimposed  they  do  not  intersect  each 
other  and  the  graphical  subtraction  cannot  be  performed.  The  Fj^rtoft 
technique  for  evaluating  the  LaPlacian  cannot  be  applied  ordinarily 
unless  large  differentiation  distances  are  used, 

7.  Possible  improvements 

In  many  applications  in  meteorology  it  is  the  distribution  of 
absolute  vorticity  about  the  vertical  which  is  desired  rather  than  the 
distribution  of  the  relative  vorticity.  This  requires  that  the  Coriolis 
parameter  be  added  to  the  relative  vorticity.  The  addition  is  frequently 


17 


omitted  with  the  justification  that  the  gradient  of  the  Coriolis  para¬ 
meter  is  relatively  small.  There  are  undoubtedly  occasions  in  which 
it  is  desirable  to  avoid  this  additional  error,  and  the  reason  for 
failing  to  avoid  it  is  the  inconvenience  incurred  with  existing  pro¬ 
cedures.  When  using  the  computer  described  here,  it  would  be  simple 
to  use  a  double  gang  switch  for  the  latitude  control,  and  have  the 
second  gang  operate  a  selective  voltage  divider  which  could  then  apply 
an  additional  potential  to  the  meter,  proportional  to  the  value  of 
the  Coriolis  parameter.  Since  the  addition  would  be  done  electrically, 
one  could  obtain  the  absolute  vorticity  or  the  relative  vorticity  at 
will,  with  no  further  effort  on  the  part  of  the  operator. 

Undoubtedly,  there  are  other  arrangements  and  designs  which  would 
be  as  effective  or  more  so  than  the  present  one.  For  example,  it  is 
possible  to  use,  instead  of  four  measuring  arms,  a  system  of  three, 
passing  through  the  vertices  of  an  equilateral  triangle.  This  would 
avoid  the  difficulty  which  is  experienced  with  the  present  apparatus 
in  the  vicinity  of  sharply  curved  contours.  The  electrical  analog 
might  then  be  more  complicated;  however,  this  question  has  not  been 
looked  into. 

Once  it  is  realized  that  any  desired  combination  of  space  deriva¬ 
tives  of  a  mapped  field  can  be  obtained  by  mechanical  measurement  of  the 
distances  between  isopleths  of  the  field  combined  with  electrical  analog 
computation,  a  whole  host  of  possibilities  suggest  themselves.  For  ex¬ 
ample ,  in  evaluating  vorticities  it  ought  to  be  possible  to  take  into 
account  the  fact  that  the  geostrophic  vorticity  is  an  overestimate  in 
troughs  and  an  underestimate  in  ridges.  The  electrical  analog  could 
be  so  constructed  as  to  include  this  automatically  in  the  result. 
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8.  Conclusions 

The  electrical  analog  computer  described  here  has  been  tested,  as 
discussed,  and  was  used  for  a  period  of  two  weeks  in  the  current  weather 
analysis  program  conducted  by  graduate  students  at  the  Florida  State 
University.  The  results  of  this  experience  indicate  that  one  can  evalu¬ 
ate  the  geostrophic  vorticity  rapidly  and  accurately  by  this  technique. 
There  is  reason  to  believe  that  the  process  of  measurement  is  suffi¬ 
ciently  simple  that  sub-professional  persons  could  carry  it  out  suc¬ 
cessfully.  With  a  little  experience,  the  procedure  should  be  faster 
than  the  previously  used  methods  without  loss  of  accuracy. 
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